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Abstract 

It  has  been  demonstrated  that  it  is  possible  to  have  a  better  understanding  of  the  mass-transfer  phenomena  occurring  at  PEMFC  electrodes  by 
distinguishing  the  various  diffusive  regimes  taking  place  inside  the  porous  layer  close  to  the  electrodes  themselves.  In  each  regime  the  interactions 
between  diffusive  and  forced  flows  have  been  expressed  in  terms  of  Peclet  numbers  and  the  overall  diffusive  resistances  have  been  expressed  in 
terms  of  Sherwood  numbers. 

Now  the  comparison  of  traditional  and  non-traditional  geometrical  arrangements  can  be  more  fully  discussed  and  the  geometrical  optimisation 
of  the  cell  can  be  more  clearly  determined  by  explicitly  considering  the  role  of  the  geometrical  arrangement  of  the  cell  channels  and  the  related 
pressure  field. 

Both  interdigitated  and  serpentine  cells  can  be  operated  in  such  a  way  that  high  Sherwood  numbers  and,  then,  a  high  limit  current  are  attained. 
Both  types  of  cells  are  penalised  by  higher  head  losses  than  traditional  cells,  but  these  appear  to  be  much  greater  in  the  serpentine  arrangement. 

The  ultimate  goal  of  reaching  high,  almost  uniform  Sherwood  numbers  and  low  head  losses  is  still  problematic.  A  partially  interdigitated 
configuration  might  be  a  step  in  the  right  direction. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

In  the  first  part  of  this  work  [1]  attention  was  focused  on 
the  gas-phase  mass  transfer  occurring  at  PEMFC  electrodes. 
The  various  possible  diffusive  regimes  were  determined  for  each 
regime,  the  interaction  between  diffusive  and  forced  flows  was 
expressed  in  terms  of  the  Peclet  number  and  then,  the  overall 
diffusive  resistance  was  expressed  in  terms  of  the  Sherwood 
numbers  as  a  function  of  the  Peclet  numbers. 

The  results  of  this  approach  have  been  used  in  the  discussion 
of  a  number  of  experimental  findings  relating  to  traditional  [2,3], 
interdigitated  [4-7]  and  serpentine  [8]  cells.  A  better  under¬ 
standing  of  the  role  of  the  gas  mass  transfer  has  been  shown 
to  be  effective  in  explaining  some  important  differences  in  the 
performance  of  these  cells,  especially  in  terms  of  limit  current. 
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Now,  in  the  second  part  of  this  work  the  comparison  of  inter¬ 
digitated  and  serpentine  geometries  will  be  further  developed  by 
explicitly  taking  into  account  the  pressure  differences  between 
contiguous  channels  and  their  dependence  on  the  head  losses 
through  the  channels.  These  results  are  currently  available  as 
the  output  of  detailed  cell  models  [7,9];  here,  a  simplified,  ana¬ 
lytical  solution  is  proposed,  by  means  of  which  the  differences 
between  different  geometrical  arrangements  can  be  immediately 
appreciated. 

Both  interdigitated  and  serpentine  cells  can  be  operated  at 
high  Sherwood  numbers  corresponding  to  high  limit  currents 
[1,9],  but  both  are  correspondingly  penalised  by  higher  head 
losses  than  traditional  cells,  this  being  much  greater  in  the  case 
of  the  serpentine  arrangement.  A  better  understanding  of  these 
features  will  allow  us  to  more  effectively  plan  the  geometrical 
optimisation  of  the  cell. 

However,  the  ultimate  goal  of  reaching  high  and  uniform 
Sherwood  numbers  in  all  parts  of  the  cell  plane  and  low  head 
losses  is  still  problematic,  and  even  interdigitated  cells  do  not 
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Nomenclature 

A 

relative  permeability — def.  (4) 

b 

width  of  a  channel  (m) 

B 

def.  (19) 

c 

characteristic  velocity  of  Eq.  (42)  (ms-1) 

Cl,  c2, 

C3,  C\  integration  constants — def.  (10) 

c5,  c6, 

C7,  Cg  integration  constants — def.  (23) 

d 

length  of  the  diffusive  layer  (m) 

df 

diameter  of  the  holes  interconnecting  the  channels 
(m) 

D 

effective  diffusivity  in  the  porous  medium 
(m2  s-1) 

Dm 

molecular  diffusivity  in  the  continuous  gas  phase 
(m2  s-1) 

/c>/f?/p 

/s  shape  factors 

h 

thickness  of  the  diffusive  layer  (m) 

k 

permeability  of  the  porous  medium  (m2) 

L 

length  of  the  channel  (m) 

M 

def.  (44) 

nc 

number  of  channels  in  a  cell 

N 

mass  flux  (kgm-2  s-1) 

P 

dimensionless  pressure 

P 

absolute  pressure  (kgm-1  s-2) 

P  r 

reference  pressure — def.  (4)  (kgm-1  s-2) 

Pex 

transversal  Peclet  number 

Pey 

axial  Peclet  number 

Sh 

Sherwood  number 

u 

velocity  in  the  channel  (ms-1) 

uv 

inlet  (maximum)  u  velocity  (ms-1) 

V 

velocity  in  the  porous  medium  (ms-1) 

vr 

reference  velocity  for  v — def.  (4)  (ms-1) 

X 

spatial  coordinate,  perpendicular  to  the  electrode 
(m) 

y 

spatial  coordinate,  parallel  to  the  electrode  (m) 

z 

dimensionless  spatial  coordinate  along  the  chan¬ 
nel 

Z 

spatial  coordinate,  parallel  to  the  channel  (m) 

Greek  letters 

A 

dimensionless  head  loss — def.  (19) 

<P 

interdigitation  degree 

viscosity  of  the  gas  (kgm-1  s-1) 

P 

density  of  the  gas  (kg  m-3) 

Superscripts  and  subscripts 

c 

channel 

f 

interconnecting  hole 

L 

end  of  the  channel 

P 

permeable  boundary 

r 

reference  value 

s 

porous  medium 

0 

beginning  of  the  channel 

0, 1, 2, 

3  four  successive  channels 

/ 

refers  to  serpentine  arrangement 

seem  to  be  able  to  supply  a  completely  satisfactory  answer  to 
this  problem  at  present.  The  partially  interdigitated  configuration 
proposed  in  the  final  section  of  this  work  might  be  a  step  in  the 
right  direction. 

Reference  schemes  for  the  geometrical  arrangements  for  an 
interdigitated  (Fig.  la)  and  a  serpentine  (Fig.  lb)  configuration, 
respectively,  considered  are  reported  in  Fig.  1. 

2.  The  flow  field  for  interdigitated  geometry 

As  discussed  in  the  first  part  of  this  work  [1],  a  good  estima¬ 
tion  of  the  velocity  of  migration  in  the  porous  medium  close  to 
the  electrode  is  an  important  and  necessary  premise  for  a  good 
estimation  of  the  Sherwood  number  and,  then,  the  limit  current. 
The  well-known  integrated  one-dimensional  permeability  equa¬ 
tion 

fMP\  -  P2) 

v=  - - -  (1) 

lid 

is  probably  accurate  enough  to  evaluate  the  velocity  of  migra¬ 
tion  when  (P\  —  P2)  <  P\  and  a  substantially  one-dimensional 
migration  path  (h^d;  see  Fig.  1,  Part  1  of  this  work  [1]),  but 
it  requires  a  knowledge  of  the  driving  force  (Pi  —  P 2)  and  this, 


Fig.  1.  Reference  geometrical  and  flow  schemes  for  an  electrode  element:  (a) 
an  interdigitated  arrangement  and  (b)  a  serpentine  arrangement. 
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in  turn,  depends  on  the  head  losses  inside  the  channels  and  the 
channel  mass  balance. 

Eq.  (1)  and  the  consequent  evaluation  of  the  local  Sherwood 
number  can  be  included  in  a  detailed  simulation  tool,  so  that  the 
effects  of  the  gas  mass  transfer  towards  the  electrodes  is  more 
correctly  taken  into  account  and  a  better  prediction  of  the  cell’s 
behaviour  at  high  currents  can  be  obtained.  Some  examples  of 
the  results  obtainable  using  this  method  have  been  reported  in 
the  first  part  of  this  work.  Here,  simplified  analytical  solutions 
for  the  flow  and  pressure  fields  are  presented,  so  that  the  dif¬ 
ferences  between  the  various  geometrical  arrangements  can  be 
fully  demonstrated. 

With  reference  to  a  pair  of  interdigitated  channels,  1  and  2, 
under  the  simplified  condition  of  constant  total  flow  rate,  the 
mass  and  momentum  balances  can  be  written 


dPi  _  fc32u\/j  dP2  _  fc32u2/J 

dZ  “  ti2  ’  dZ  “  b2 

du\  fvvh  du2  fvvh 
~dZ  ~  W'  ~dZ  ~+  b2 


(2) 

(3) 


where  each  channel  is  assumed  to  be  in  a  laminar  regime  and 
the  migration  flow  rate  lost  by  the  first  channel  is  acquired  by 
the  second. 

The  constant  volumetric  flow  rate  assumption  (in  synthe¬ 
sis  u\  +  U2  =  constant)  needs  some  further  explanation.  Strictly 
speaking,  it  should  correspond  to  the  assumption  of  null  cur¬ 
rent,  because  the  electrochemical  reaction  necessarily  implies  a 
flow  rate  variation,  a  decrease  at  the  anode  and  an  increase  at 
the  cathode.  Moreover,  any  temperature  distribution  in  the  cell 
is  associated  with  further  volumetric  flow  rate  changes.  Bearing 
in  mind  these  considerations,  the  reference  to  a  mean,  constant, 
volumetric  flow  rate  could  seem  an  oversimplification.  On  the 
other  hand,  the  pressure  drops  to  be  considered  are,  primarily, 
those  taking  place  at  the  cathode,  where  the  molar  flow  rate  vari¬ 
ations  of  the  supplied  air  are  of  the  order  of  a  4-6%,  so  that  the 
assumption  of  a  mean  flow  rate  value  implies  quite  acceptable 
uncertainties  (less  than  3%).  Errors  on  the  anodic  side  could 
be  significantly  greater  (namely  30-40%),  but  the  problem  of 
the  anodic  pressure  drops  is  significantly  less  important.  How¬ 
ever,  if  one  is  interested  in  a  more  exact  evaluation  of  the  flow 
field  on  both  sides,  taking  effective  flow  rate  changes  as  well 
as  temperature  distributions  into  account,  detailed  cell  models 
are  available  [7,9].  The  simplified  approach  used  here  offers 
the  advantage  of  an  explicit  analytical  solution  that  is  easily 
suited  to  a  substantially  quantitative  discussion  of  different  cell 
arrangements. 

By  combining  Eqs.  (l)-(3)  and  by  defining  the  dimensionless 
parameters  and  variables 


vr  = 


uYb 2 

7p hL9 


fc32ur/jiL 
b 2 


/c/P/s64ML2 

db 4 


(4) 

(5) 


the  following  system  of  differential  equations  is  obtained: 


d2(pi-p2 )  a2( 

- -  =  A  (pi-  p2), 


d2(pi  +  pj) 

dz2 


=  0 


For  an  interdigitated  geometry  the  boundary  conditions  to  be 
added  to  the  equations  in  (6)  are 


z  =  0;  u\  =  uY\  ii2  =  0, 

dp2  _  d(p\  -  P2)  _ 

dz  ’  dz 


z  =  i; 

u  i  =  t 

j;  u  2  =  uY , 

dp2  

-i, 

<N 

1 

dz 

dz 

The  solution 


Pi  ~  Pi  =  C i  exp(Az)  +  C 2  exp(-Az), 

spi  +  P2  =  C3z  +  C4  (9) 


with 


1  +  exp(-A) 
exp(A)  —  exp(— A) 

1  +  exp(A) 
exp(A)  —  exp(— A) 


(10) 


demonstrates  that  the  migration  velocities  v,  which  are  propor¬ 
tional  to  the  pressure  difference, 


v 

vY 


l)m  —  VY 


(11) 


are  symmetrical  to  the  middle  of  the  cell  (z=  1/2),  while  the 
velocities  along  the  channels,  which  are  proportional  to  the  pres¬ 
sure  derivatives,  are  anti- symmetrical.  The  overall  head  losses 
for  the  pair  of  interdigitated  channels  is 

Ci[l  +  exp(A)]  +  C2[l  +  exp(-A)]  +  1 
P 10  -  P2L  =  - 2 - 

and  the  fraction  of  the  entire  flow  rate  that  migrates  through  the 
porous  medium  from  the  first  channel  to  the  second  (by  pass 
fraction)  corresponds,  obviously,  to  unity  (vm/vY  =  1). 

The  main  parameter  determining  the  flow  field  configura¬ 
tion  is  A,  which  is  a  function  of  the  permeability  of  the  porous 
medium  and  the  geometry  of  the  system.  The  effects  of  parame¬ 
ter  A  on  the  flow  fields  and  the  pressure  field  in  an  interdigitated 
arrangement  are  illustrated  in  Fig.  2. 

By  varying  A  the  following  limit  cases  can  be  considered. 
When  the  permeability  is  much  higher  than  the  head  losses  of 
the  channels 
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Fig.  2.  The  distribution  of  the  flow  rate  between  two  adjacent  interdigitated 
channels  along  the  axial  coordinate  for  various  values  of  parameter  A:  (a)  the 
velocity  inside  the  first  channel  (u/ur)  and  (b)  the  migration  velocity  between 
channels  (v/vT).  (c)  The  trend  of  the  global  head  losses  (p\o  —  P2if)  as  a  function 
of  parameter  A. 


•  A»  1 

mi  1  +  exp(— Az)  —  exp(— A  +  Az) 

-  - 

X  2  (13) 

j  [exp(-Az)  +  exp(— A  +  Az)] 

1 

P 10  -  P2L  ~  ~  (14) 

half  of  the  entire  flow  rate  is  almost  immediately  transferred 
from  one  channel  to  the  other,  then  each  channel  conveys  half 
the  flow  rate,  with  overall  head  losses  almost  equal  to  a  half 
of  the  ones  pertinent  to  a  single  channel  with  a  velocity  uY  (or 
equal  to  the  ones  pertinent  to  a  single  channel  with  a  velocity 
uY/ 2).  The  remaining  half  flow  rate  is  transferred  to  the  second 
channel  immediately  before  the  exit  point. 


v  (A 

_  I  _ 

Fr  ~  \2 


Obviously,  a  pair  of  interdigitated  channels  does  not  oper¬ 
ate  in  this  way.  The  correct  operating  conditions  are,  on  the 
contrary,  nearer  to  the  second  limit  case,  where  the  perme¬ 
ability  of  the  porous  medium  is  low. 

•  A«1 


u  i  v 

—  ~  1  —  z,  —  ~  I 

UY  VY 


P 10  -  p  2L  ^ 


(15) 

(16) 


and  the  overall  head  losses  are  controlled  by  this  low  perme¬ 
ability  and  are  much  greater  than  the  ones  of  a  single  channel. 
In  this  case,  the  velocity  of  migration  v  is  almost  uniform 
along  the  entire  channel  length  (z  axis). 

As  an  effective  interdigitation  should  require  A  <  1 ,  the 
global  head  losses  of  the  cell  are  correspondently  high  (see 
Fig.  2c).  A  good  compromise  can  be  reached  by  choosing 
values  of  A  of  the  order  of  3,  so  that  a  good  distribution  of 
the  interdigitation  flow  is  achieved  in  correspondence  to  the 
global  head  losses  that  are  only  a  little  greater  than  those  of 
the  traditional  arrangement. 


3.  The  flow  field  for  serpentine  geometry 


The  serpentine  geometry,  also  in  its  simplest,  single  channel 
version,  is  a  little  more  difficult  to  analyse  than  the  interdigitated 
one.  The  study  of  two  contiguous  channels  (1  and  2),  now  with 
countercurrent  velocities,  needs  to  consider  another  two,  exter¬ 
nal,  channels  (0  and  3),  which  are  also  involved  in  an  exchange 
of  mass  with  the  preceding.  On  the  other  hand,  the  simplification 
of  the  repeatable  trend  of  the  head  losses  along  the  serpentine 
channel  can  be  introduced 


Po(z)  =  P2(z) +  2A,  p3(z)  =  Pi(z)  -  2A 


(17) 


and,  then  a  formulation  strictly  similar  to  (6)  can  be  attained 


d2(pi  -  pi-  A) 

dz2 

d2(p  i  +  pi) 
s  &  =° 


B2(Pi  ~  Pi-  A), 


(18) 


B2  =  2 A2,  2 A  =  POL  -  PlL  =  P2L  -  P3L  =  P 10  -  P20 

(19) 


On  the  other  hand,  the  correct  set  of  boundary  conditions  is 
now 


z  =  0;  U2  =  ~u  i,  p\-p2  =  2A, 

d(pi  +  pi)  = 
dz 


(20) 


Z  =  1;  «2  =  -Ml,  Pl~P2  =  0, 


d(pi  +  pi) 
dz 


=  0 
(21) 
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so  that  the  solution 

pi  -  P2  =  C5  exp (Bz)  +  C(5  exp (—Bz)  +  A, 

Pi  +  P2  =  C7z  +  C8  (22) 


c  _  A[ 1  +  exp(— ff)] 

5  exp(£)  —  exp(— B) 

c  _  A[  1  +  exp(£)]  (23) 

6  exp(£)  —  exp(— B) 

C7  =  0 


demonstrates  that  this  time  u  is  symmetrical  and  v  anti- 
symmetrical  to  the  middle  of  the  cell: 

U]_  _dp]_  _  _f  1\  d(gi  -  /72) 

wr  dz  V  2  /  dz 

_  /  A£\  [l+exp(-£)]  exp(£z)+[l  +  exp(£)]  exp(-£z) 
\  2  /  exp(Z?)  —  exp(— B) 

(24) 


v 

vY 


B< 


(Pi  ~  P2), 


Vm  <  Vi 


(25) 


The  by-pass  fraction,  corresponding  to  the  mean  value  of  v, 
is  now 


v 


m 


B- 


Vr 


(Pi  ~  P2)m  = 


B< 


A 


(26) 


and  the  sum  of  (24)  and  (26)  corresponds  to  unity,  that  is  it  rep¬ 
resents  the  whole  reference  flow  rate  through  the  cell.  This  last 
condition  makes  it  possible  to  evaluate  the  head  loss  parameter 


A  = 


B < 


+ 


B 


exp  (B)  +  exp(— B)  +  2 
exp(Z?)  —  exp(— B) 


(27) 


The  effects  of  parameter  A  on  the  flow  fields  and  the  pressure 
field  in  a  serpentine  arrangement  are  illustrated  in  Fig.  3. 

With  regard  to  the  limit  cases,  when  the  permeability  is  high 

•  £»  1 


A 


B2  ’ 


vm  ~  v, 


(28) 


(C)  A 


Fig.  3.  The  distribution  of  the  flow  rate  between  two  adjacent  serpentine  chan¬ 
nels  along  the  axial  coordinate  for  various  values  of  parameter  A  (A2  =  B2/ 2):  (a) 
the  velocity  inside  the  first  channel  ( u/ur )  and  (b)  the  migration  velocity  between 
channels  (v/vf).  (c)  The  trend  of  the  global  head  losses  (pio  ~P2o)  as  a  function 
of  parameter  A. 


•  5«1 


A  ss 


1. 


% 

vr 


« 1 


U 1 
UY 


1, 


z)  <  1 


(30) 

(31) 


—  ~  \  )  [exp(— 5z)  +  exp(— 5  +  Bz)\ 

wr  \B  J  ^29) 

—  ^  1  +  exp(-5z)  -  exp(-Z?  +  Bz) 
vY 

almost  the  entire  flow  rate  is  transferred  through  the  porous 
layer  (the  by-pass  fraction  is  only  a  little  less  than  unity),  while 
the  velocity  along  most  of  the  channels  and  the  corresponding 
flow  rate  fraction  are  very  low. 

This  first  limit  case  may  seem  interesting,  but  it  is  not  very 
realistic:  even  in  the  serpentine  geometry  the  most  likely  limit 
case  is  the  one  corresponding  to  low  permeability. 


which  controls  the  overall  head  losses.  In  this  case,  the  veloc¬ 
ity  of  migration  v  is  very  low  and  strongly  asymmetrical, 
while  the  head  losses  of  the  channel  are  only  slightly  affected 
by  the  very  low  by-pass  fraction. 

In  a  serpentine  cell  the  single  channel  head  losses  are 
higher  than  in  the  traditional  arrangement  because  the  veloc¬ 
ity  is  higher;  the  decrease  corresponding  to  the  increase  in 
the  by-pass  fraction  (see  Fig.  3c)  is  not  decisive;  more¬ 
over,  the  global  head  losses  depend  on  the  number  of  pas¬ 
sages,  so  that  increasing  the  cell  size  quickly  makes  them 
unacceptable. 
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4.  The  order  of  magnitude  of  the  parameters 

We  can  assume  unit  values  for  the  various  shape  factors  and 
we  have  referred  to  literature  values  for  the  geometrical  and 
physical  parameters: 

h  —  2  x  1 0  2  cm,  L  =  10  cm,  d  =  b  =  1 0  1  cm 

&=10-8cm2,  fi  =  10-4  gem-1  s-1, 
sAPr  =  3  x  103  Pa  =  3  x  104  g  cm-1  s-2 

Dm  =  0.2  cm2  s  ,  - =0.1 

Du 

and  we  have  also  defined  [1] 


cases.  Assuming  low  permeability  regimes  (A  <<C  1;  £  <<C  1),  the 
interdigitated  case  is  characterised  by 


/c64wr/xL 

A2/?2 


while  the  serpentine  is  characterised  by 


Un 


B2vt 


B2u[b 2 
fphL  ’ 


\AP\  = 


fc32u'r/iL' 

b2 


Under  the  same  cell  flow  rate  and  the  same  number  of  chan¬ 
nels  for  the  cell, 

,  uTnc 


ur  = 


L  =  Ln( 


the  comparison  of  the  by-pass  fractions 


the  axial  Peclet  number:  Pey  =  dv/D 
the  transversal  Peclet  number:  Pex  =  h2v/dD 
and  the  Sherwood  number  Sh  as  the  ratio  between  the  actual 
mean  flux  of  the  electrochemical  reagent  and  a  reference  flux 
[1]. 

Then,  by  making  these  assumptions  the  transition  from  the 
diffusive  to  the  forced  regime  can  be  situated  around  Pey  =  25  and 
the  transition  from  the  forced  flow  rate  to  the  forced  pellicular 
regime  around  Pey  =  60.  Moreover,  the  above  set  of  data  leads 
to 

A2  =  0.13,  A  =  0.35,  £  =  0.5 


(Un)serp  B2n c  A2nc  (32) 

Unterd  8  4 

even  when  not  considering  the  highly  asymmetrical  distribution 
of  the  serpentine,  strongly  favours  the  interdigitated  geometry, 
unless  the  number  of  channels  is  very  high.  But  in  such  a  case 
the  comparison  of  the  head  losses 


|A£|Serp  A2n2 

A  P  |  interd  4 


(33) 


should  demonstrate  that  the  overall  serpentine  head  losses  are 
much  greater  than  the  already  considerable  ones  that  characterise 
an  interdigitated  cell. 


so  that  such  approximations  as  exp(A)  ~  1  +  A  are  correct  to 
within  a  few  percent  and,  then,  a  sure  reference  can  be  made 
to  the  limit  cases  of  low  permeability.  However,  it  should  be 
kept  in  mind  that  a  relatively  little  change  in  the  geometrical  or 
physical  parameters  is  also  sufficient  to  make  this  simplification 
less  accurate  and  raise  the  values  of  B  to  the  order  of  unity. 

For  the  interdigitated  geometry,  with  A  =  0.35,  the  velocity  of 
migration  and  the  corresponding  Peclet  numbers  are 

n  =  30cms_1,  £<^  =  150,  Pex  =  6,  Sh  =  0.9 

So,  the  migration  through  the  porous  medium  takes  place 
under  a  forced  regime  and,  prevalently,  under  a  forced  pellicular 
regime.  Under  such  conditions  the  Sherwood  number  is  of  the 
order  of  unity. 

For  the  serpentine  geometry,  with  £  =  0.5,  the  velocities  of 
migration  are  undoubtedly  lower;  a  mean  value  is  about 

v  =  2  cm  s_1 ,  £^  =  10,  Pex  =  0.4,  Sh  =  0.3 

so  that  the  migration  still  takes  place  in  a  substantially  diffusive 
regime,  with  Sherwood  numbers  of  the  order  of  0.3. 

A  significant  increase  in  the  velocity  of  migration  and  the 
Sherwood  number  can  only  occur  when  the  head  losses  are 
strongly  increased,  as  is  clearly  demonstrated  by  the  results 
reported  by  Zhukovsky  and  Pozio  [8]. 

Another  way  of  looking  at  the  question,  so  that  the  differ¬ 
ent  flow  conditions  of  a  serpentine  channel  and  an  interdigitated 
configuration  can  more  clearly  be  taken  into  account,  is  to  con¬ 
sider  the  by-pass  fractions  and  the  overall  head  losses  in  the  two 


5.  The  partial  interdigitated  geometry 

In  the  preceding  section,  we  demonstrated  how  the  serpentine 
geometry  could  contribute  to  lowering  the  diffusive  resistances 
only  if  heavy  overall  head  losses  are  allowed.  On  the  other  hand, 
Zhukovsky  and  Pozio  have  underlined  how  a  significant  lower¬ 
ing  of  the  diffusive  resistances  can  also  be  attained  by  using 
geometries  that  are  characterised  by  rather  moderate  migration 
velocities. 

The  interdigitated  geometries  also  have  greater  head  losses 
than  the  traditional  geometries  but,  on  the  other  hand,  they  permit 
migration  velocities  that  are  often  greater  than  what  is  needed. 
For  instance,  if  the  migration  velocity  of  the  above  example  were 
of  20  cm  s-1  instead  of  30  cm  s_1 ,  we  would  still  have  Pey  =  100, 
a  forced  flow  rate  regime  and  Sh  =  0.6,  sufficient  for  a  significant 
increase  in  the  limit  current. 

A  useful  modification  to  the  interdigitated  geometry  is  what 
could  be  referred  to  as  “partially  interdigitated”:  instead  of  feed¬ 
ing  the  entire  flow  rate  of  a  pair  of  channels  only  to  the  first 
channel,  so  that  it  transversally  migrates  to  the  second  and  exits 
entirely  from  that  one,  one  could  feed  a  fraction  (1  +  cp)/ 2  of  the 
flow  rate  to  the  first  channel  and  a  fraction  (1  —  cp)/. 2  to  the  sec¬ 
ond,  so  that  only  the  fraction  cp  (interdigitation  degree)  is  forced 
to  migrate  transversally  from  one  channel  to  the  other  and,  at 
the  exit,  the  first  channel  still  conveys  the  fraction  (1  —  <p)!2  to 
the  exit  and  the  second  the  fraction  (1  +  cp)/ 2.  The  correct  distri¬ 
bution  of  the  flow  rate  in  the  two  channels  could  be  obtained  by 
means  of  calibrated  holes  and  (low)  localised  head  losses. 
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The  differential  equations  in  (6)  are  still  valid  for  this  partial 
interdigitated  geometry,  but  the  boundary  conditions  have  to  be 
changed  in  order  to  allow  for  the  different  distribution  of  the 
flow  rate 


dpi 

d  z 


d (pi  ~  Pi) 
d  z 


dp  i 

d  z 


d(p\  -  pi) 

d  z 


1  +  (p 


(34) 


(35) 


A  solution  such  as  (9)  is  still  obtained,  with  the  integration 
constants  C\  and  C2  multiplied  by  a  factor  (p  over  the  values  in 
(10)  and,  again,  C3  =  —  1.  The  velocities  of  migration  through 
the  porous  layer,  which  are  still  symmetrical,  are  reduced  by  the 
factor  cp 

j  (p\  -  Pi),  Vm  =  (pVr  (36) 


just  as  the  overall  head  losses  of  the  two  channels  are  similarly 
reduced 


P 10  -  P2L  = 


\  [1  +  exp(A)][l  +  exp(-A)]  l 
A)  exp(A)  —  exp(— A)  2 


(37) 


The  by-pass  fraction  obviously  corresponds  to  the  interdigitation 
degree  (vm/vY  =  cp).  In  particular,  the  low  permeability  limit 
case  now  predicts 


•  A«  1 

u\  1  +  cp  —  2  cpz 

-  - 

uY  2 

2cp 

P 10  “  P2L  ~  J2 


V 

—  ~  (p 

VY 


(38) 

(39) 


so  that  the  overall  head  losses  can  be  significantly  reduced 
without  hampering  the  limit  current  too  much. 


The  partially  interdigitated  geometry  could  present  a  further 
advantage.  As  has  already  been  said,  partial  interdigitation  can 
be  realised  by  means  of  calibrated  holes  that  connect  a  pair  of 
channels  at  their  beginning  and  at  their  end  so  that  the  flow 
rate  distribution  between  the  two  depends  on  the  localised  head 
losses  of  the  holes 


1  2 


<P) 


Recalling  Eq.  (1)  and  making  reference  to  low  permeabil¬ 
ity  (A  <  1),  so  that  (Pi  -  ^2)0  =  (Pi  -  P2)l  =  Pi  -  P2  =  AP  and 
v(z )  =  constant,  we  have 


v  =  cp 


uYb2 

hd 


MAP 

pid 


and  the  following  equation  in  cp  can  be  written: 


<P  = 


(1  -  (p2)ur 


c  = 


2  pddf 


c  ffpklb 2 

The  interdigitation  degree 


<P  = 


—  1  +  [1  +  4(«r/c)2]1/2 
2  uY/c 


(41) 


(42) 


(43) 


is  then  an  increasing  function  of  the  velocity  uY.  In  turn,  the 
migration  velocity 


v  =  McpuY,  M  = 


b* 


hd 


(44) 


shows  a  more  than  linear  increase  with  uY ,  that  is  the  whole 
flow  rate  of  the  two  channels.  So,  when  the  desired  current, 
and  then  the  flow  rate,  is  increased,  the  utilisation  degree  being 
equal,  the  partially  interdigitated  cell  responds  with  an  increase 
in  the  interdigitation  degree  and  a  more  effective  increase  in 
the  migration  velocity,  the  Sherwood  number  and  the  limit  cur¬ 
rent,  at  the  expense  of  a  lower  increase  in  the  overall  head 
losses. 
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